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VO Abstract. We present a new systematic method to construct Abelian functions on Jacobian 

varieties of plane, algebraic curves. The main tool used is a symmetric generalisation of 
the bilinear operator defined in the work of Baker and Hirota. We give explicit formulae 
for the multiple applications of the operators, use them to define infinite sequences of 
Abelian functions of a prescribed pole structure and deduce the key properties of these 
functions. We apply the theory on the two canonical curves of genus three, presenting new 
explicit examples of vector space bases of Abelian functions. These reveal previously unseen 
similarities between the theories of functions associated to curves of the same genus. 
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^ 1 Introduction 

O 

The overall aim of this paper is the presentation of a systematic method of construction of 

• • Abelian functions on Jacobian varieties of plane, algebraic curves. Apart from its purely theo- 

.£h retical interest this is important from the point of view of at least two applications. Firstly, 

the relations between Abelian functions of given pole order, that is, the expansion of given 

functions in terms of the basis, are partial differential equations of the kind arising in integrable 

systems theory. In particular, members of the KP hierarchy have been shown to arise in this 

way, for functions associated with hyperelliptic curves of general genus and also for various non- 

hyperelliptic curves. Consequently, the Abelian functions constitute one class of solution to such 

partial differential equations. This connection is explicit even in the early work of Baker (see 

later references) . Secondly, Abelian functions with shifted argument expanded in terms of the 

basis of functions in un-shifted arguments give addition laws generalizing the classic addition 

law of points on the cubic curve. This expresses the group structure of the Jacobian. It is also 

important in the theory of discrete integrable systems. 

The need for a new systematic method of construction is motivated in Section 1.2 by the so 

called basis problem for Abelian functions. The solution we present is based on our definition 

of generalised Baker-Hirota operators, following the work of Baker on Abelian functions and 

the work of Hirota in soliton theory. Hence we continue the introduction in Section 1.1 with 

a discussion of the Hirota derivative and how it motivates the general operators which follow. 
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1.1 Hirota derivatives and operators 

The Hirota derivative is originally a device belonging to the direct method of solution to soliton 
equations, [27]. A simple paradigm is provided by the Korteweg-de Vries equation, 

Ut + u xxx + 6uu x = 0. 

Substitution for u(x, t) by 2(lnF(x, t)) xx and integration with respect to x (neglecting arbitrary 
functions of t) yields the bilinear form 

D x (D t + D 3 x )F-F = 0, 
where the Hirota derivative is defined as 

D n x D?F ■ G = 9£d?F(x + x',t + t')G(x -x',t- OL=o,f =0" 

Thus, for example, 

F) X F • G = F X G — FG X , 

D x D t F ■ G = F xt G - F x G t - F t G x + FG xt 

etc. The iV-soliton solutions arise via an ansatz of the form F(x, t) = f(e ei , e 6 * 2 , . . . , e &N ) where 
the 0i are linear in x and t. 

Of course, the Hirota "derivative" is not actually a derivative since it does not act on products 
but on pairs as signalled by the • above. However, if we consider a product on function pairs of 
the kind (f,g)(h, k) = (fh,gk) and define T> x etc. as Hirota operators that map from function 
pairs to function pairs, 

Vx(f,g) = (d x f,g) - (f,d x g), 

then we do have a Leibniz derivation formula: 

V x ((/, g)(h, k)) = V x (f, g)(h, k) + (/, g)V x (h, k). 

Because further, this Hirota operator is bilinear over addition and scalar multiplication, it is 
sensible to treat the function pairs as tensor products and the product as the standard multi- 
plication on tensor products: 

(f®g)(h«S)k) = fh®gk. 

The original Hirota derivative is then a composition of T> x with a total symmetrization opera- 
tor, S, mapping from the tensor to the differential polynomial algebra: 



F®F 


Vx 


F®F 


s 


F, 


f®g 


T> x 


fx®g-f®g x 


s 


fxg - fg x 



Here F denotes a 6-algebra of appropriate ^-valued functions, for some field t. 

Such a generalisation has been the subject of previous work by one of the authors, applied 
to invariant theory [2], the theory of Pade approximants [1] and to generalised p- functions [5]. 

It is also interesting to compare with algebraic notions of derivations over tensor products of 
algebras introduced in [31]. It is not clear to the present authors whether T> is either an example 
or a generalization of this approach but it suggests an avenue for further investigation 

In each of these applications, there is an equivariance property of the Hirota operators ap- 
propriate to the case. Suppose we have a group, Q acting on the space of independent variables 
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represented by the list x. We assume there is an induced representation of Q on some set of 
functions of x: 

f{G- 1 x)=f G (x). 

In some applications Q is a Lie group and we are interested in a Lie algebra representation. 
Then a Hirota operator T>, (we drop the subscript), should have the desirable property: 

F&F 4 F(g>F 

F®F 4 F(g>F 

In the case of the Abelian functions associated to a curve, just such an action is obtained. Given 
a curve of genus g we define the associated Abelian functions as those meromorphic functions of 
g complex variables u G C 9 , which are periodic with respect to the period lattice of the curve. All 
such functions may be expressed algebraically using the cr-function of the curve. A translation 
across lattice points, u t- ► u + ■£ induces a multiplicative factor: o(u + £) = h(u,£)a(u). The 
generalised Hirota operators we introduce here will generate Abelian functions out of tensor 
products of these a-functions. 

Our generalisation of the Hirota operator picks up from earlier work. In [24], trilinear Hirota 
derivatives are presented in the context of certain integrable equations in trilinear form. They 
are characterised by a simple equivariance property. Further work can be found in [25, 26]. 

In [2] a more algebraic, independent treatment is developed for Hirota type operators acting 
on arbitrary tensor products of finite-dimensional sl2(C)-modules. In the infinite dimensional 
limit the 5(2 action becomes a Heisenberg action and the Hirota operators become the familiar 
Hirota derivatives (when acting on degree two tensor products). This interplay has also been 
exploited in a study of differential normal forms, [34] . 

In the present paper we deal with arbitrary tensor products of a differential algebra over 
some field, £, which may as well be C given the analytic context. The differential algebra is 
generated by a Kleinian u-function in variables u\, ... ,u g and its partial derivatives. There is 
a Heisenberg action on the derivatives with respect to each individual U{ and an s^C) action 
on the whole set of U{. For example the set of first order derivatives (d Ul a, d U2 a, . . . , d u a) is 
a basis for a g dimensional s[2(C)-module. 

1.2 The basis problem for Abelian functions 

The main motivation for the generalised Hirota operators presented in this paper is to solve the 
so called basis problem for Abelian functions. The simplest Abelian functions are the elliptic 
functions, associated with elliptic curves, which of course have genus one. Elliptic functions 
have been the subject of much study since their discovery and have been extensively used to 
enumerate solutions of non-linear wave equations. They occur in many physical applications; 
traditionally the arc-length of the lemniscate and the dynamics of spherical pendulums, [32], 
but also in cryptography, [36], and soliton solutions to the KdV equation, [14]. 

Recent times have seen a revival of interest in the theory of their generalisations. These 
Abelian functions are also beginning to find a wide range of applications. For example, they 
give further solutions to the KdV equation along with solutions to other integrable equations 
from the KP-hierarchy (see for example [9, 11, 21]). They have also been used to describe 
geodesic motions in certain space-time metrics, [23]. 

Much progress has been made through the realisation of Abelian functions as generalisa- 
tions of the Weierstrass p-function, following the approach of Klein and Baker for hyperelliptic 
functions. In the genus g case the original Weierstrass p-function is generalised to a family of 
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functions denoted pij(u) for i,j £ (l,...,g). It is the case that py = pji and so there are 
\g{g + 1) functions in total. We have the integrability property that 

d d 

So we can denote the derivatives of the functions using additional indices and adopt the notation 
of writing the indices in ascending numerical order. The 3-index p- functions generalise the 
derivative of the Weierstrass function, p' and so on. 

We may classify the Abelian functions according to their pole structure. The m-index p- 
functions will have poles of order m for example. We can define a vector space of Abelian 
functions with poles of at most a given order and so the question of building bases for the vector 
spaces naturally arises. The Riemann-Roch theorem for Abelian varieties gives the dimension of 
the vector space as m 9 where m is the maximal pole order. In the genus one case each successive 
vector space will have dimension one higher, and the new basis function can be played by the 
next derivative of the p-function. However, for g > 1 there are not enough p-functions to build 
the bases. Indeed, as both g and m increase so does the deficit of functions. This is the basis 
problem. 

Solutions to the problem have been found in a number of specific cases. In genus 2 the 
introduction of one further function with poles of order 3 is sufficient to allow differentiation 
to generate all subsequent bases. This new function is the difference P11P22 — Pi2> m which the 
poles of order 4 in each term cancel. In [16, 19] the basis problem was solved in various cases, by 
a group including one of the present authors, through the introduction of analogous functions; 
polynomials in p with coefficients chosen to cancel poles. However for curves with g > 2 
hyperelliptic and g > 3 non-hyperelliptic the bases cannot be finitely generated by differentiation 
and this approach does not generalise easily to give an infinite number of functions. 

Solving the basis problem makes it easier to work with the functions in a wide range of 
applications. For example, in [22] a reduction for the Benney equation was constructed using 
the Abelian functions associated with a genus six curve and having a basis for the functions 
limited the search for the explicit form of the reduction. The bases can also be used to derive 
more theoretical results, such as the differential equations satisfied by the functions. These 
can be derived from a set which generalise the classic differential equation for the Weierstrass 
p- function, 

(p'(u)) = Ap{uf - g 2 p{u) - g z . 

The bases can also be used to construct various addition formulae for the functions. For exam- 
ple, there are generalisations of the following well known addition formula for the Weierstrass 
functions, 

a(u + v)a(u - v) . , . . 

a(u) ji a{v) z 

Each generalisation has the same ratio of <r-functions expressed as a polynomial in an appropriate 
basis of Abelian functions. For curves with an extra cyclic symmetry in their coordinates 
there are other classes of addition formulae. For example, suppose the curve is invariant under 
[C J ] : (x,y) —> (x,C 3 y) where £ is some primitive mth root of unity. Then we can define 
a corresponding action on the variables u and find an expression in an appropriate basis of 
Abelian functions for 



n 
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where u*- 1 ', . . . , w m ' are different sets of variables. The existence of such formulae helped moti- 
vate the new Abelian functions constructed in Section 6. For examples of using bases to calculate 
differential equations and addition formulae see [9, 16, 18, 19, 21]. 

While the p-functions may be defined algebraically using the curve, an alternative definition 
using the <T-function of the curve can be easier to work with. Recall that Weierstrass introduced 
an auxiliary function, cr(u), in his theory which satisfies 



■2 



d 

The Weierstrass a- function is usually defined as an infinite product over the periods, but it can 
also be expressed using the first Jacobi ^-function, multiplied by a constant and exponential 
factor. This definition generalises naturally to give higher genus a-functions as multivariate 
functions defined using the Riemann ^-function Once the generalised cr-function is defined we 
may then consider generalised p-functions in analogy to the previous equation: 

Pij (u) = -g^ j lo g [a(u)]. 

Then, using the notation above, we see this is equivalent to the following definition in Hirota 
operators: 



fci = ( 2~~2 ) s oV i o'D j (a® a) 



where T>i indicates that the differentiation in the Hirota operator is with respect to u%. This led 
to the definition of so called Q- functions, where the operator is applied not twice but an even 
number of times, generating further Abelian functions with poles of order 2. See [9, 18, 20, 21] 
for examples of how such functions have been used to complete pole 2 bases. In this paper we 
provide the first proof that such functions are always Abelian. 

These Q-functions led to the derivation of the generalised Hirota operators which are intro- 
duced in this paper. They allowed for the definition of more general class of Abelian functions 
which we call "/^-functions. Infinitely many such functions can be defined with a given pole 
order, making them useful for the general basis problem. 

The paper is divided into two parts. The first deals with the theory of the operators and the 
second the use of this theory in the construction of Abelian functions. While it is this application 
which is of primary interest in this paper, it is anticipated that the operators may have other 
uses, particularly in the world of integrable systems, and so we present their theory separately 
where possible. Section 2 reformulates the original Hirota theory into a form that more easily 
generalises, and proves some explicit results for the functions obtained by repeated application 
of the operators. Section 3 then defines the generalised operators and derives some of their 
properties. Sections 5 and 6 apply the operator theory to define infinite classes of Abelian 
functions, first with poles of order 2 and then poles of arbitrary order. Finally Section 7 uses 
these functions to construct new bases of Abelian functions for the two canonical curves of genus 
three. 

2 The Baker— Hirota operator 

We start by fixing some notation. Throughout the paper we are dealing with functions of g 
variables. We drop the variables from our notation unless required. So for example, 

/ = f(u) = f(u 1 ,U 2 ,...,Ug). 
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When working with Abelian functions, g will be the genus of some underlying algebraic curve. 
Until we come to the work in Section 7, none of the calculations in this document will depend 
on the particular curve or genus. In this section and the next the reader may assume we are 
just working in a C-algebra of entire functions and their tensor products (^) c . 

Later we use subscripts on a function to indicate the application of operators, but for now we 
shall use subscripts to indicate differentiation with respect to one of the variables. For example, 



d_ 



fi = ^r/^i'^, ...,%). 



Superscripts in square brackets indicate a label, with other superscripts indicating the object 
raised to a power as normal. 

Definition 2.1. We let T>i denote the Baker-Hirota operator. This acts on the tensor product 
of two functions as 

This operator is the key component in the famous Hirota derivative, discussed in the intro- 
duction. We refer to it as a Baker-Hirota operator since it was first used by Baker in [8]. 
We note that Baker-Hirota operators commute with each other, 

v i ov j (fin®fW)=v j ov i {fm®f®). 

Definition 2.2. We let S denote the C-linear symmetrization operator which will send a tensor 
product to standard multiplication, 

S : / ® g — > fg = gf. 

We are usually interested in applying the Baker-Hirota operator to a tensor product of the 
same object. We note that the symmetrization operator would annihilate the resulting sum after 
the application of a single operator, 

SoVi(f<8)f) = S(fi<8>f-f<8> fi) = fif - ffi = 0. 

However, when two operators are applied, there is a non-trivial result left after symmetrizing, 

SoV i0 Vj(f ® /) = S^fij ® / - fi ® fi) - (fj ® fi - f ® fa)) = 2(ffij - fifj). 

We will now consider the effect of multiple applications of a Baker-Hirota operator. We use 
In = {h,i2, ■ ■ ■ , in} for the set of indices of the operators, assumed distinct for the moment. We 
introduce the multiple sum notation, 

n n n 

E -££-£■ 

l<jl<J2<---<jm<n jl=lJ2=jl+l jm=jm-l + l 

Note that if there are no values of an index j in the appropriate range then that sum does not 
run. 

Lemma 2.3. For n > 1 we have 

V il o...oV in {f®g) = JT (-l) m Y. hn\{i h ,i j2 ,.,i im } ® 9ini j2 ...i jm ■ (2-1) 

m=0 \l<ji<J2<— <j m <n J 

When m = the entire inner sum is interpreted as a single term with the set {ij 1 ,ij 2 , ■ ■ ■ ,ij m } 
being the empty set, 0. A function with an empty set of subscripts is just the function itself 
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We give the details of an inductive proof in Appendix A. 

The lemma allows us to observe the following result, detailing when the symmetrization 
operator annihilates the result of Baker-Hirota operators applied to the same function. 

Corollary 2.4. If n is odd then S o V ix o • • • o V in (/ <g> /) = 0. 

Proof. We will show that when n is odd, symmetrizing in equation (2.1) causes all terms to 
cancel. For example, the first term, fj n <g> /, and the last term, (—1)™/ ® fj n will cancel with 
each other. 

First note that if n is odd then the outer sum in (2.1) will have an even number of terms. 
We may consider these terms pairwise, when m is mi = [§J — k and 1712 = [§] + k, for some 
integer k. So each pair consists of the sums 



\ I / j •> 1 n\\1j 1 ^j 2 >—'hm 1 i J hl' 

l<jl<J2<--<jm<n 



'J2>""' 'Jrni 



and 



(-1) 



m-2 



£ 



JIr l \{ij 1 , 



l<jl<J2<-~<jm<n 



'32''""' 'imj 



}®4-i. 



'J2'""'*3r»j| 



Now, n is odd so the floor, [J and ceiling |~] of § will be of opposite relative parity, and hence 
so will ttt-1 and rri2 as they differ from the floor and ceiling by the same integer k. 
Now since 



n = 
we have 

n — m\ 



n 




[W 


— 


+ 


— 


V'l- 




2 



/ 


n 




vn~ 


\ / 


n 


, \ 


rni 




— 


+ 


— 


- 


— 


-k = 


— 


V 


L2_ 




2 


) V 


L2_ 


/ 


2 



+ k = 1712, 



and similarly n — ni2 = mi. Hence the partition of indices on either side of the tensor in the 
two terms are complementary. When we symmetrize they will give the same sum of terms, and 
so the differing signs outside the sums will cause them all to cancel. ■ 

Similar calculations will verify that the object is non-zero when n is even. 
We now introduce a multiplication on tensor products, which simply creates a tensor product 
of the same length by multiplying corresponding entries. For example, 

(a <8> b) (c <g> d) = (ac) (bd) . 
The operation is clearly symmetric; 

(o(8)6)(c(8)d) = (ac) 8) (bd) = (c<g> d)(a (g> b). 
Lemma 2.5. The Baker-Hirota operator satisfies a Leibniz rule. That is 

A((a0 6)(c(»d)) = (a®b)Vi(c®d) + (c®d)Vi(a®b). (2.2) 

Proof. By direct calculation we see that both sides of (2.2) are 

ctjC ®bd + aci bd — ac<S> bid — ac® bdi. ■ 
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Corollary 2.6. The Baker-Hirota operators satisfy a general Leibniz rule, 

n 

V h o-.-o V in ({a®b){c®d)) =J2J2v 7Tl (a(g)b)V n2 (c(g)d). 

k=0nGU 

Here U is the set of ordered pairs of disjoint partitions, it, of the set of indices I n into subsets, 
7Ti and 7T2 of length n — k and k respectively. The symbol T>- Ki represents the concatenation of 
Baker-Hirota operators with indices the entries in TX{. 

Proof. Recall that the standard product rule (fg)' = f'g + fg' leads by induction to the general 
Leibniz rule, 



U9) in) =jz{t)f {k) 9 { 

l — n \ / 



where (Tj are the binomial coefficients. This leads identically from Lemma 2.5 to the following 
result for a single Baker-Hirota operator, 

v^({a ®b){ C ® d)) = Y, uWfc ® b ) v r k ( c ® d )> 

where the superscripts on the operator mean it is applied multiple times. If instead we apply 
different Baker Hirota operators then we must consider all the different possibilities for k ap- 
plications of an operator in each step. This leads to the definition of LT as in the corollary. The 
number of different possibilities is the number of ways of choosing k from n, and so still (?) . ■ 

3 Generalised Baker— Hirota operators 

In this section we consider a generalisation of the Baker-Hirota operators which maintains many 
of the properties from the previous section. Such operators will be classified with an order m, 
and will be multi-linear, acting on a tensor product of length m. We first introduce an operator 
to differentiate within such tensor products. 

Definition 3.1. Define of to be an operator which differentiates the jth entry in a tensor 
product with respect to the variable Ui, i.e. 

of (/W ® ■ ■ ■ (8) /b'- 1 ! /0 8) f l3+1] ® • • • ® / H ) 
= f [1] <g> • • • <g> f [j ~ 1] 8) fl j] 8) f [j+1] <g> • • • ® f [m] . 
We note that these operators commute with each other, 



ab'ilflb'a] _ aih] a\jl] 
For brevity we will often express these tensor products using a product notation, such as 

m 



fc=i 

We now define the generalised operators. We note that these operators are similar to the 
(symmetric operators) proposed in [24]. However, the authors there considered only the trilinear 
cases (m = 3) in detail, and actually preferred to use a different non-symmetric basis that could 
be described as products of bilinear operators. Symmetric multilinear Hirota maps of the sort 
defined below were discussed in [3] in the context of algebraic covariant and invariant theory. 
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Definition 3.2. Define / H i to be the mth order generalised Baker-Hirota operator. This acts 
on the tensor product of m functions as 

m m m 

yH. ®/w->x;c i - i fl{ , ' i (8)/ w . 

fc=l 3=1 k=l 

where C is a primitive mth root of unity and m > 2. 

We can take £ to be any primitive mth root of unity, but we may assume without loss of 
generality that 

C = exp . 

\ m J 

Recall that 

i + C + C 2 + --- + C m_1 = o. (3.1) 

We note that the 2nd order generalised Baker-Hirota operator is exactly the original Baker- 
Hirota operator from Definition 2.1, 

%f ] = Vi. 

For brevity we may write only "^-operator unless we mean to distinguish between the cases. 

Since the df -operators commute with each other, it also clear that they commute with any 
^-operator, 



From this we can conclude that the %-operators of the same order commute with each other: 



V'i=l / h=l 

We extend the symmetrization operator from Definition 2.2 to apply it to tensor products of 
arbitrary length, 

m m 

fe=l fc=l 

Again, we are usually interested in applying the generalised operators to a tensor product of the 
same object. As before, the symmetrization operator would annihilate the resulting sum after 
the application of a single operator: 

(m \ / m m \ m 

fc=i / \i=i fc=i / j=i 

where the final equality follows from equation (3.1). In this general case we find the sym- 
metrization operator will also annihilate the result after two applications, so long as m > 2. For 
example, consider applying the 3rd order operator twice: 

3 

,[3] 



W x l ® / = (°fh ®f®f + ( 1 f®fn®f + ( 2 f®f®fh, 



k=l 
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3 

«!? ° nl S ® / = C%,i 2 ® / ® / + CVn ® /i 2 €5 / + C 2 /n ® / ® fi 2 
fc=l 

+ CVt a ® /n ® / + C 2 / ® /iaA ® / + C 3 / ® /n ® /i a 

+ C 2 /i 2 ® / ® fh + C 3 / ® /*, ® /u + C 4 / ® / ® JW (3.2) 

So on symmetrizing we have 

3 

s °hJ 31 owg ] (g)/ = / ni2 / 2 (C° + C 2 + C 4 ) + /u/fa/CC 1 + C 2 + C 1 + c 3 + c 2 + c 3 ). 

fc=l 
Then since £ 3 = 1 we have 

3 



r[3] „ a/t 3 ] 

fe=l 



5 o ftg J o «gJ (g) / = f llt2 f{l + c 1 + C 2 ) + 2/ n / J2 /(l + c 1 + c 2 ) = o. 



To see what happens in general we start by considering the effect of multiple applications of 
the operators. 

3.1 Multiple applications of the generalised operators 

We are not able to use the multiple sum notation of Lemma 2.3. Instead we introduce some 
new notation, which may seem a little cumbersome at first, but is useful for proving the later 
results. This notation will be used throughout the rest of the paper. 

Once again let /„ = {i\, 12, ■ ■ ■ , i n } be a set of indices. Let P™ be a set whose elements p are 
sequences constructed from the partitions of the number n of length at most m. If a partition 
has length less than m then we extend it to length m by adding zeros. We denote by p\, . . . , p m 
the parts of the partition in decreasing order and then we define p to be the sequence of parts 
[pi, . . . ,p m ]. Next we let ^(p) be the set of permutations, ip, of the sequence p. We denote by 
■01, • • • , ipm the entries in a permuted sequence ip. We define a function z(ip) which returns an 
integer according to 

m 

z(0) = X>-l)W-)O- (3-3) 

fc=i 

Finally we let Tl(ip) be the set of all disjoint partitions tt of the set of indices I n into m subsets 
of lengths given by the sequence ip. We denote the resulting subsets of such a tt by m, ■ ■ ■ , ir m - 

Lemma 3.3. For n > 1 we have 



in 



*T ° • • • ° wT ® flk] = £ £ c * w £ ® 'S ■ ( 3 - 4 ) 

fc=i peP™ i>e9(p) 7ren(v>) fc=i 

Proof. We do not need to give an inductive proof of Lemma 3.3 as it may be determined directly 
from a careful analysis of the resulting differential operator acting on the tensor product, 

m I n I m 

fc=i \fc=i y=i 

The polynomial in <9-operators would be fully symmetric if we had £ = 1. We organise the terms 
as described in the lemma so that we collect together those with the same coefficient in £ and 
we observe that this coefficient is then given using the function in equation (3.3). ■ 
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We note that the size of the set n(/0) is the number of ways of splitting n elements into m 
groups of sizes ipi, . . . ,ip m . This is by definition a multinomial coefficient. Further, since ip is 
just a permutation of p, we have that 



|n(V0l = |n(p)| ' " ' 



P!,...,p m J Pl\---p 



!...o„! 



rn ■ 



We now give an example of the construction to demonstrate the notation. Consider the m = 3, 
n = 2 case that we looked at earlier. Here I n = {41,12} an d so P$ are the partitions of 2 with 
length at most 3. There are two partitions, and we extend them each to length 3 to give 

p [1] = [2,0,0] and p [2] = [1, 1, 0]. 
First we consider p^ 1 '. There are three permutations of [2,0,0] and so three entries in \l/(pW); 

^ [1] = [2,0,0], V [2] = [0,2,0], and ip [3] = [0,0, 2]. 

We see that 

,(#]) = (0)(2) + (l)(0) + (2)(0) = 0, 

and similarly z(t/j^- 2 ') = 2 and z{il)™) = 4. For each ip G ^ there is only one it G n(-0) since there 
is only one possible split of the indices into subsets of these lengths. We have 

n(#]) = {[{h,i 2 }, {},{}]}, nfo&M) = {[{}, {il> hh { }]} 
and n(^ [31 ) = {[{},{}, {kM}}- 

Hence p'- 1 ' contributes three terms to the sum; 

cvli ® / [2] ® / [3] + c 2 / [i] ® 41 ® / [3] + c 4 / [i] ® / [2] ® /5 



'_> 



Now we consider p™. There are also three permutations of [1, 1, 0] and so again three entries in 

#1 = [1,1,0], V [2] = [1,0,1], and V [3] = [0,1,1]. 

We calculate z(ip^') = 1, z(^™) = 2, z(V>™) = 3. In each of these case there are two possible 
partitions of the indices into sets of these lengths. For the first we have 

n(V [11 ) = {[M,M,{}], [{i 2 },M,{}]}. 

This contributes the terms 

c 1 (/ff®/£ I ®/ M +/ff®^ I ®/ M ) 

to the sum. Similarly we have Il(^ 2 ') = { [{h}, {}, {22}], [{^2}, {}, {*i}] } contributing 

c 2 (ti ] ®f l2] ®fff+fl>f l2] ®fl?) 

and II^) = { [{}, M, {t 2 }], [{}, M, {ti}] } contributing 

c 3 (/ [11 ^/|>/1 2 31 + / [1 ^/1>/1 31 ). 

So we see we have found the nine terms expected from our calculation in equation (3.2). 
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3.2 Application on the tensor product of a function with itself 

We will now restrict to the case where the functions in the tensor product are taken to be the 
same, and consider the effect of the symmetrization operator. In this section we make use of 
the theory of symmetric functions. We gather together the necessary definitions and theory in 
Appendix B. 

Lemma 3.4. Let X = [X±, . . . , X m ] be given by X^ = C fc ~ 1 . Then for n > 1 we have 

m / m \ 

So«H ... oW Mg )/= £ Mp{x) £ (HuX (3.5) 

fe=i peP™ 7ren(p) \fc=i / 

where M are the monomial symmetric functions {see Definition B.l). 

Proof. We consider equation (3.4) from Lemma 3.3. Under the symmetrization operator S the 
tensor product becomes a regular product, so 




7ren(V>)A:=l 



n* 



Now if all the functions p> = f then this sum over H(ip) is the same for all ip G ^(p)- So 
we may instead denote it as a sum over 11 (p) and factor this entire object outside of the sum 
over ^(p). We cannot discard the permutations ip entirely though, as the term f"' still varies, 
even under symmetrization. We have 

m I / m 

so*Ho...o«H(g)/=x; e n/- 

fe=i p€P™ \7ren(p) \fe=i 
We further analyse the sum in the final factor. We have 

m 

e ^ )= e n^- 1 )^- 




If we introduce the dummy variables Xk = ( then we have that 

m 

E < ZW = E \[xt h = M p {X u ..^X m \ 

from Definition B.l. Hence the final factor is identified as M p {X) as described in the lemma. ■ 

Denote the power sum symmetric functions by pk (see Definition B.2). When evaluated on X 
as in Lemma 3.4 we find (see Lemma B.3) that 

{m ifm|fc, 
otherwise. 

By considering expressions in the power sum basis we can determine when multiple applications 
of the ^-operators map to zero. 

Corollary 3.5. Ifm\n then 

m 

SoH^ ] o...oH t ™ ] (g)f = 0. 
fc=i 
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Proof. Every monomial symmetric function is homogeneous in degree. For those appearing in 
Lemma 3.4, the degree is always n. We may think of the degree as a weight. Since each power 
sum pk has degree k, the expression for M p as a polynomial in pi, . . . ,p m will only involve terms 
where the subscripts of the power sums add to n. Hence if m \ n then the expression cannot 
contain any terms with p m as the only variable. However, it is only these terms that do not 
evaluate to zero with the choice of variables X^ = Q k ~ 1 by Lemma B.3. 

So if m\n then every term in the right hand side of equation (3.5) is zero. ■ 

Corollary 3.6. If m\n then S o W™" o • • • o 9i\ ®fcLi / is n °t identically zero. 

Proof. Consider the expression in Lemma 3.4. One of the terms in equation (3.5) will cor- 
respond the partition on n with p = [n, 0, . . . ,0]. But for such a partition M p = p n and by 
Lemma B.3 we know that this power sum is not zero for m\n. So one of the terms in the 
expression is a non-zero constant times a monomial in derivatives of / which does not occur in 
any of the other terms. Hence the expression cannot be identically zero. ■ 

We can go further and give a formula for the value of the M p in the case where n\m. To 
do this we use a procedure of Doubilet from [15] to express the monomial symmetric functions 
in power sums with coefficients given by Mobius functions, /^(vr). We give the full details in 
Appendix B with the key formula stated in equation (B.2). 

Corollary 3.7. The augmented monomial symmetric functions on X = [1, £, . . . , £ m_1 ] are 

M P =Y1 tin)™', 
Trerr 

where £ is the number of subsets in ir and II' contains only those set partitions tt S n such that 
all parts ofv(p,ir) are multiples of m. 

Proof. From Lemma B.3 we know that when evaluating equation (B.2) on the variables 
Xk = £ , only those terms such that all parts of v(p,ir) are multiples of m will be non- 
zero. In these cases each pi is m and so each Mobius function has a coefficient of m to the 
power of the length of v(p, 7r). ■ 

The monomial symmetric functions involved in Lemma 3.4 will evaluate to the integers for 
the augmented functions from the corollary above, divided by the constants from equation (B.l). 
However, we observe that when m\n then Sow™ o • • ■oW ' ®a!Li / a ^ so nas integer coefficients 
as a polynomial in the derivatives of /. In fact, these integer coefficients all have a common 
factor of m. 

3.3 Leibniz rule 

We recall the multiplication operation we defined for tensor products, which multiplies corre- 
sponding entries in tensor products of the same length. We use it again here to show that the 
%-operators all satisfy Leibniz properties. 

Lemma 3.8. An mth order %-operator satisfies a product rule. That is 




(3.6) 
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Proof. Using the definitions we have 

m m 

LHS(3.6) = ^C J '- 1 #(8)/ [ S [fcl , 

3=1 
m I / m \ / m \ / m \ / m 



[J] 





RHS(3.6) = £ <T X ® /W a?"' g)»W + 0/] dj«g)/W 

i=i V \k=i / V fe=i / \/c=i / V *:=i / 

Comparing the jth term in each sum shows them to be equal, (and hence that the operator d, 
satisfies its own product rule). 

Corollary 3.9. The %- operators satisfy a general Leibniz rule, 

n / m 

e=o vren \fc=i / 

Here IT is again the set of disjoint partitions, tv, of the set of indices I n into two subsets, 7Ti 
and 7T2, of length n — £ and I respectively. The symbol 7~L™ represents the concatenation of 
%-operators of the same order with indices the entries in 7Tj. 

Proof. This follows the proof of the original Leibniz rule and the proof of Corollary 2.6. ■ 

We note that as in Corollary 2.6, the size of II is still the binomial coefficient ("). 

4 Abelian functions associated with algebraic curves 

The first generalisation of the elliptic functions to hyperelliptic functions was pioneered by Klein 
and Baker as described in Baker's classic texts [6, 7] . They defined multivariate p- functions as 
discussed in Section 1. Hence these generalised functions are sometimes called Kleinian. The 
more general definitions presented on this section are primarily from the work of Buchstaber, 
Enolskii and Leykin in [11] and the further developments in [17] and [33]. 

The periodicity conditions of all the higher genus functions are defined using the period 
lattice of an underlying algebraic curve. We use the higher genus a- function of the curve to 
define them. We give explicit definitions relating to the class of (n, s)-curves, defined below. 
This class is special as the curves all have a branch point at infinity. The restriction to this class 
is so that we can make use of the well developed theory of series expansions for the functions 
associated to such curves to check the linear independence of functions. 

However, we note that one can define Abelian functions associated to curves outside this 
class. One of the authors has pioneered a so called equivariant approach to the functions in [4], 
where the curve and functions all transform under an sl2 action. Another recent development 
is in [28] where the authors define a-functions associated to Riemann surfaces without making 
any realisation of the surface as an algebraic equation. The definitions and results in Sections 5 
and 6 for Abelian functions defined using %-operators should be applicable in these alternative 
theories. It is only in the calculation of the explicit bases in Section 7, that we make use of the 
(n, s)-theory. 

Definition 4.1. For two coprime integers (n, s) with s > n we define an (n, s)-curve, denoted C, 
as an algebraic curve defined by f(x, y) = 0, where 

f(x,y) = y n + qi^y 11 - 1 + q 2 (x)y n - 2 + ■■■ + q n - 1 (x)y - q n (x). (4.1) 
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Here x, y are complex variables and qj(x) are polynomials in x of degree (at most) [js/n\ . We 
define a simple subclass of the curves by setting qj(x) = for < j < n — 1. These simpler 
curves are defined using curve parameters Xj as 

f{x, y) = y n - (x s + Xs-ix 8 ' 1 + • • • + Aix + A ) (4.2) 

and are called cyclic (n,s)-curves. In the literature the word 'cyclic' is sometimes replaced by 
'strictly' or 'purely n-gonal'. Note that we also have curve parameters for the general (n,s)- 
curves, contained within the qj(x). 

Restricting to the cyclic classes results in much easier computations, but it not usually ne- 
cessary for theoretical reasons. However, the cyclic curves do possess extra symmetry which 
manifests itself in the fact that the associated functions satisfy a wider set of addition formulae. 

If the curve C is a non-singular model of X then the genus is given by g = ^{n — l)(s — 1). 
The associated functions will be multivariate with g variables, u = (ui, . . . , u g ). For example, 
a (2,3)-curve would have genus one and is indeed an elliptic curve. The associated Weierstrass a 
and p-functions depend upon a single complex variable u. 

4.1 Defining the functions 

We first describe how the period lattice associated to the curve may be constructed. We start 
by choosing a basis for the space of differential forms of the first kind; the differential 1-forms 
which are holomorphic on the curve, C. There is a standard procedure to construct this basis 
for an (n, s)-curve (see for example [33]). In general the basis is given by hidx/f y , i = 1, . . . ,g 
where the hi are monomials in (x, y) whose structure may be predicted by the Weierstrass gap 
sequence, although other normalisations are sometimes used. 

We next choose a symplectic basis in Ui(C, Z) of cycles (closed paths) upon the compact 
Riemann surface defined by C. We denote these by {a±, . . . ,a g , )3i, . . . ,(3 g }- We ensure the 
cycles have intersection numbers 

aiCtj = 0, PiPj = 0, ai(3j = 5ij 

The choice of these cycles is not unique, but the functions will be independent of the choice. 

We introduce dr as a basis of differentials of the second kind. These are meromorphic 
differentials on C which have their only pole at oo. This basis is usually derived alongside the 
fundamental differential of the second kind. Rather than repeat the full details here we refer 
the reader to [6] for the general theory and [10] which gives a detailed example construction. 

We can now define the standard period matrices associated to the curve as 

u/ = I d> dui 

\Ja k / k,f=l,...,g 




rf = I ® dr e 



u = | 


f duA 


■;9 


7]" = ( 


f drA 

^Jfa / k,e=i,.. 


;9 



'a k / k,£=l,...,g 

and define the period lattice A formed from uj',uj" by 

A= {oj'm + oj"n, m,n£Z 9 }. 

This is a lattice in the space C 9 . The Jacobian variety of C is presented by C 9 /A, and is denoted 
by J. We define k as the modulo A map, 

k : C 9 ->■ J. 



16 M. England and C. Athorne 

For k = 1, 2, . . . define 2lj, the Abel map from the k-th. symmetric product Sym (C) of C to J 
by 

2l t : Sym fc (C) -»• J, 

(Pi,...,P fc ) M- ( / du+-+ du) (mod A), 



OO JOO 



where the p are again points upon C. Denote the image of the k-th Abel map by W*- k > and 
define the k-th standard theta subset (sometimes referred to as the k-th strata) by 

e [k] = w [k] U [-l]W [k] , 

where [—1] means that 

[-1] («!,...,%) = (-ui,...,-u g ). 

We are considering functions that are periodic with respect to the lattice A. 

Definition 4.2. Let 9Jt(u) be a meromorphic function of u £ C 9 . Then SDt is a standard Abelian 
function associated with C if it has poles only along k -1 (0^ -1 ]) and satisfies, for all £ £ A, 

£Dt(u + £) = 9Jt(w). (4.3) 

All the Abelian functions in this paper are defined using the c-function of the curve. While 
not an Abelian function itself, it does satisfy a quasi-periodicity property. Let 8 = oj'8' + uj" 8" 
be the Riemann constant with base point oo. Then [S] is the theta characteristic representing 
the Riemann constant for the curve C with respect to the base point oo and generators {ocj, Pj} 
of #i(C,Z). (See for example [11, pp. 23-24].) 

Definition 4.3. The Kleinian a-function associated to a general (n, s)-curve is defined using 
a multivariate ^-function with characteristic S as 



a(u) = cexp (luri'iu'y^eiS}^')- 1 ^ \ {J)- l u") = cexp {\un\J) 



- l u T ) 



x Y, ex P 2-Ki{\{m + 6') T {J)- l J'{m + 8') + {m + 8') T {{J)- l u T + 8")} 

The constant c is dependent upon the curve parameters and the basis of cycles and is fixed 
later, following Lemma 4.8. 

We now summarise the key properties of the a- function. See [11] or [33] for the construction 
of the cr-function to satisfy these properties. For any point u £ C 9 we denote by u' and u" the 
vectors in ~R 9 such that 

U = U! U + U U . 

Therefore a point £ £ A is written as 

e = j a + uj"£" £ a, £', i" £ u>. 

For u, v £ C 9 and £ £ A, define L(u, v) and x(^) as follows: 

L(u,v) = u T (n'v' + rf'v"), 

i\Tsii i e"\T X ' i l(e'\To'n 



X (£) = exp [2ri{{£') T 5" - (£") T 5' + \{£') T £"}]. 
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Lemma 4.4. Consider the a-function associated to an (n,s)-curve. 

• It is an entire function on C 9 . 

• It has zeros of order one along the set k~ 1 (@^ 9 ~ 1 '), commonly called the Q-divisor. Further, 
we have cr(u) ^ outside the set. 

• For all u G C 9 , £ G A the function has the quasi-periodicity property: 

£ 



a(u + £) = x(-^)exp 



L[ u + 2 



a(u). (4.4) 



We let L = L(-u + |) and h = \ eL denote the extra factor introduced. So a(u+£) = ha(u). 
• It has definite parity given by 

a{-u) = (-l)^™ 2 - 1 ^ 2 - 1 ^). 

Proof. The function is clearly entire from the definition, while the quasi-periodicity and the 
zeros are classical results (see [6]), that are fundamental to the definition of the function. They 
both follow from the properties of the multivariate ^-function. The parity property is given by 
Proposition 4(iv) in [33]. ■ 

We can now define p-functions using an analogy of the relation between p and a in the 
elliptic case. 

Definition 4.5. Define m-index Kleinian p-functions for m > 2 by 

do o 

Ph,i 2 ,...,im{u) = -q T, o lOg O-(ti) , 

du h du i2 du im 
where i\ < ■ ■ ■ < i m G {1, . . . , g}. 

The m-index p-functions are meromorphic with poles of order m when cr(u) = 0. We can 
check that they satisfy equation (4.3) and hence they are Abelian by substituting equation (4.4) 
directly into their definition. The m-index p-functions have definite parity with respect to the 
change of variables u — > [— l]u. This is independent of the underlying curve, with the functions 
odd if m is odd and even if m is even. The ordering of the indices is irrelevant and so for 
simplicity we always order in ascending value. 

The new Abelian functions introduce in Sections 5 and 6 are defined using the operators 
applied to tensor products of functions. We define an operator for checking their periodicity. 

Definition 4.6. Given the lattice A, define a corresponding lattice operator by 

Ai : u \— > u + £, 

for £ some arbitrary point in the lattice. So Ai(o~) = a(u + £) = ha for example. In a slight 
abuse of notation we also use A^ for the operator acting on a tensor product, 

n n 

A^(g)/ H =(g)A£/ [il . 

i=\ i=\ 



in 



Lemma 4.7. The operator A^ commutes with S, V>i and H.^ . 
The proofs are all simple by direct calculation. 
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4.2 Weights and expansions 

For a given (n, s)-curve we can define a set of weights, denoted by wt and often referred to as 
the Sato weights. We start by setting wt(x) = —re, wt(y) = —s and then choose the weights 
of the curve parameters to be such that the curve equation is homogeneous. We see that for 
cyclic curves this imposes wt(Aj) = — n(s — j). The Abel map 2ti gives an embedding of the 
curve C upon which we can define £ = x" as the local parameter at the origin, 2li(oo). We 
can then express y and the basis of differentials using £ and integrate to give series expansions 
for u. We can check the weights of u from these expansions and see that they are prescribed by 
the Weierstrass gap sequence. 

By considering Definition 4.5 we see that the weight of the p-functions is the negative of 
the sum of the weights of the variables indicated by the indices, irrespective of what the weight 
of o~(u) is. We note that curves of the same genus will, notationally, have the same p- functions, 
but may exhibit different behavior as indicated by the different weights of the variables and 
functions. We will discuss the weight of the cr-function below. All other functions discussed are 
constructed from a or p-functions and their weights follow accordingly. We can show that all 
the equations in the theory are homogeneous in these weights, with a more detailed discussion 
of this available, for example, in [21]. 

We can construct a series expansion of the cr-function about the origin, as described below. 

Lemma 4.8. The Taylor series expansion of o~(u) about the origin may be written as 

oo 

a(u) = KSW nyS (u) + Y,C k (u). 

k=0 

Here K is a constant, SW n)S the Schur- Weierstrass polynomial generated by (n,s) and each C^ 
a finite, polynomial composed of products of monomials in u of weight k multiplied by monomials 
in the curve parameters of weight —(wt(a) — k). 

Proof. We refer the reader to [33] for a proof of the relationship between the cr-function and 
the Schur- Weierstrass polynomials and note that this was first discussed in [12]. We see that 
the remainder of the expansion must depend on the curve parameters and split it up into the 
different Cj, according to the weight split on terms between the parameters and u. We can see 
that each C/% is finite since the number of possible terms with the prescribed weight properties 
is finite. In fact, by considering the possible weights of the curve coefficients we see that the 
index k in the sum will, in the cyclic case, actually increase in multiples of n. ■ 

The Schur- Weierstrass polynomials are Schur polynomials generated by a Weierstrass par- 
tition, derived in turn from the Weierstrass gap sequence for (re, s). See [33] and [12] for more 
details on these polynomials. This connection with the Schur- Weierstrass polynomials allows 
us to determine the weight of the cr-function as 

wt(cr) = (l/24)(n 2 -l)(s 2 -l). 

In Definition 4.3 we fix c to be the value that makes the constant K = 1 in the above lemma. 
Some other authors working in this area may use a different constant and in general these choices 
are not equivalent. However, the constant can be seen to cancel in the definition of all Abelian 
functions, leaving results between the functions independent of c. Note that this choice of c 
ensures that the Kleinian cr-function matches the Weierstrass cr-function when the (re, s)-curve 
is chosen to be the classic elliptic curve. 

The expansion can be constructed by considering each Ck in turn, identifying the possible 
terms, forming a series with unidentified coefficients and then determining the coefficients 
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by ensuring the expansions satisfy known properties of the a- function. For example, using its 
vanishing properties or the identities between the associated Abelian functions. Large expansions 
of this type were first introduced in [10], which used the generalised cr-function to construct 
explicit reductions of the Benney equations. Since then they have been an integral tool in 
the investigation of Abelian functions. Recently computational techniques based on the weight 
properties have been used to derive much larger expansions and we refer the reader to [21] 
and [20] for a more detailed discussion of the constructions. 

Such expansions are possible for the general (n, s)-curves, but the calculations involved are 
far simpler for the cyclic cases. In this paper such expansions are used only to check linear 
independence in the examples of Section 7. The number of Ck required to do this will depend 
on the weight of the functions in question. The expansion will need to go far enough to give the 
monomials in curve coefficients that may be present in any identity between functions. 

4.3 Bases of Abelian functions 

We can classify the Abelian functions according to their pole structure. We denote by T(m) the 
vector space of Abelian functions defined upon J which have poles of order at most m, occurring 
only on k _1 (0I 9_1 1), where the 9 and <7-functions have their zeros and the Abelian functions 
their poles. This is usually called the B-divisor. 

A key problem is the generation of bases for these vector spaces. Note that the dimension of 
the space T(m) is m 9 by the Riemann-Roch theorem for Abelian varieties (see for example [29]). 
The first step in constructing a basis is to include the entries in the preceding basis for T{m — 1). 
Subsequently only functions with poles of order exactly m need to be sought. The Kleinian p- 
functions are natural candidates and indeed are sufficient to solve the problem in the elliptic 
case. Here g = 1 and so only one new function is required at each stage, which can be filled by 
the repeated derivatives of the Weierstrass p-functions. 

However, if g > 1 then new classes of functions are required to complete the bases. For 
example, in the genus 2 case we find that when considering T(3) we need an additional entry 
after including all possible p-functions. The function 

A = P11P22 - P12 

is usually taken to fill this hole. The individual terms in A have poles of order 4, but when taken 
together they cancel to leave poles of order 3 as can easily be checked using Definition 4.5. We 
then proceed to consider T(4) and find that two new functions are needed after the inclusion of 
the 4-index p-functions. The two derivatives of A can play this role. Naturally the derivatives 
of the A have poles of order four but we must also ensure that they are linearly independent of 
the other functions. This can be checked trivially by noting that they are of a different weight 
to both each other and all the other elements of the basis. 

It is simple to check that when considering higher values of m, the preceding basis and its 
unique derivatives always give the required number of functions for a new basis, and that they 
are all of unique weight and so linearly independent. Hence the general basis is as described in 
Table 1, where {•} is indicating all functions of this form. 

The genus 1 and 2 cases, in which we get to a stage where new bases are calculated from the 
old ones and their derivatives, are special. They fall into the class where the theta divisor is non- 
singular and the P-module structure of such cases is discussed in [13]. For g > 2 hyperelliptic 
and g > 3 non-hyperelliptic curves this is not the case and so new methods are required to derive 
bases. 

One approach to the problem is to define new functions by matching poles in algebraic 
combinations of p-functions so that they cancel, analogous to A above. This method was 
applied successfully to derive bases for various vector spaces in [19] and [16]. However, this 
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Table 1. Tabic of bases for Abelian functions associated with a genus two curve. 



Space 


Di 


mension 


Basis 




r(o) 

r(i) 
r(2) 
r(3) 

r(4) 




1 

1 
4 
9 
16 


{1} 
{1} 

{l,pll,pl2, P22> 

{1, Pll, Pl2, P22, Pill, P112, P122, P222, A} 
{1, . . . , A, pun, . . . , ^2222, 9lA, <9 2 A} 




T(m) 




2 

m 


{basis for T(m - 1)} U {{p ix ...i m } , {d h ■ ■ ■ d im . 


- 2 A}} 



approach does not generalise easily, which led to the development of the new classes of Abelian 
functions in Sections 5 and 6. 



5 Q-functions 

In this section we study the Q-functions, defined using the original Baker-Hirota operator, 
proving in particular that they are all Abelian. The results of this section also act as a motivation 
for the more general case discussed in Section 6. 

Definition 5.1. We define the n-index Q-functions from the Kleinian cr-function as 

Qh,...,in = (-^2) {S°V h o...oV in (a®a)), i x < ■ ■ ■ < i n G {1, . . . ,<?}. (5.1) 

We note that the indices on the Q-functions do not indicate differentiation, but rather the 
application of operators. The n in the definition is just the number of operators applied, and 
is not related to the n in an (n, s)-curve. Indeed, the results of this section and the next are 
independent of the underlying curve. This definition is equivalent to the following one more 
commonly found in the literature, using the Hirota derivative from the introduction 



Qi 1 ,i2,..-,in\ u ) — r, I \9 »1 *2 



Di„a(u + v)cr(u 



v=0 



We prefer Definition 5.1 as the notation generalises in a much simpler way. 

We use the label Q to follow the notation of Baker in [8], who defined the 4-index functions. 
The general definition above was developed in [21] and [20] when 6-index functions were needed 
to explicitly solve basis problems. Baker likely chose the notation Q since the functions can be 
viewed as a generalisation of the 2-index p-functions. Comparing with Definition 4.5 we see 



Q 



1.1 



(7- 



We aim to prove that the Q-functions are Abelian. So we need to show that for £ £ A 



V»i, 



(« + €) = Q iu i 2 ,..., in (u). 



We can check the periodicity of a 2-index Q-function by simply substituting in the quasi- 
periodicity property of cr(u) to the expression in u-derivatives above. While this approach 
can be used to prove the periodicity of Q-functions for other values of n, it does not extend to 
an inductive proof due to the symmetrization operator, S, in the definition. 
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Corollary 2.4 shows us that the n-index Q-functions are zero when n is odd. A similar 
calculation shows that for n even the functions are divisible by two. Lemma 2.3 may be used 
to give the general formula for an n-index Q-function in cr-derivatives and we may show by 
induction that 

<7h,i 2 ,...,i n (u + £) = xe L £M Yl a I n \{in^-^m} L ^i L in ' ' ' L ^m ( 5 - 2 ) 

m=0 \l<ji<h<—<jm<n J 

using the fact that second derivatives of L(u) are zero. However, substituting these properties 
into the formula for the Q-functions and observing the cancelations is not trivial. 

Instead, we use the properties of the operators to prove the result using the following property 
of h, the extra factor introduced in the quasi-periodicity condition (4.4). 

Lemma 5.2. We have S o V^ o ■ ■ ■ o T>i n {h <S> h) = for all n > 1. 

Proof. We note that if n is odd then the result follows directly from Corollary 2.4. To show 
the result also holds for n even consider Lemma 2.3 with f = h, 

sov il o...ov in {h®h) = f>ir( E ^\{^ 1 ^,..^}^ 1 ^,-^)»( 5 - 3 ) 

Since all the second derivatives of L are zero we have hi = x eL Li, hij = x eL L{Lj, and by 
a simple induction, 

h iu ... jin =xe L L il ---L in . (5.4) 

Note that every term in the inner sum in equation (5.3) above contains each index in /„ exactly 
once. Hence each of those terms are equal to x 2e2L ^h ' ' ' Li n . The number of terms is the 
number of ways of choosing m from n, the binomial coefficient. Hence 

SoV il o-..oV in (h®h)= X 2 e 2L L h ■ ■ ■ L in ^ (-l) m Q . 

m=0 ^ ' 

The lemma then follows from 

0. ■ 



s>i>'"(:) 



m=0 

We now use this result to give a simple proof of the periodicity of the Q-functions. 

Theorem 5.3. The n-index Q-functions are all Abelian. 

Proof. By Definition 4.6 of the lattice operator we see that the Abelian property is equivalent 
to showing 

A<?o Q ilt ... t i n = Qi lt ...,i n . 

Using the definition of the Q-functions in equation (5.1), Lemma 4.7 which concluded that A^ 
commutes with the other operators, and the periodicity property of a, we have 



^£°Qh,...,i n = [~ 2h^a 2 ) S ° Vil °'"° Vin ° M "® ")' 



(5.5) 



22 M. England and C. Athorne 

where h = \e L '. We observe that under the multiplication operation for tensor products, 

A £ (a <8> a) = {ha <g> her) = {h <8> h)(a <g> a). 
So we may use the Leibniz rule from Corollary 2.6 to give 

n 

V n 0-.-0 V in {h®h){a®<j) = ^2^V ni {h ® h)V n2 (a ® a). (5.6) 

fc=o wen 

When k = we see there is only one n G II and that tv\ is the empty set. So the inner sum in 
equation (5.6) produces just one term, 

(h(g)h)V il o...oV in {<j®a). 

When k = 1 there are n entries in II and in each case i\\ is of length one. Hence the inner sum 
in equation (5.6) produces n terms 

V h (h ® h)V In \ h {a®a)+ V l2 (h <g> h)V In \ i2 (<r ® a) + • ■ • + V in {h ® /i)P /n \ in (o- ® a). 

Every term here involves a Baker-Hirota operator applied to (/i ® /i) and hence under sym- 
metrization each term becomes zero by Lemma 5.2. Similarly, for all k > the inner sum in 
equation (5.6) produces terms which all involve k operators applied to (h ® h). Only the term 
from k = will not vanish under symmetrization. Hence we have 

S o V h o • • • o V in ((h ® h)(a <g> o-)) =li 2 So%o...o £> in (o- <g> <r) 

and so equation (5.5) becomes 

A^oQij j„ = ( ~ 2h 2 a 2 ) h2 S° v ii °---oV in (o®o) =Qi 1 ,..., in . ■ 

We have already noted that the 2-index Q-function is exactly the function pij. It is actually 
possible to express any Q-function as a polynomial in the p- functions. In [18] the authors showed 
that 

Qijkt = pijke — Zpijpke — 2pikpj£ — ZpuPjk, 
and similarly in [21] that 

(°iijklmn — Pijklmn ^ \PijPklmn "T PikPjlmn "T pilpjkmn t Pimpjkln ~r pinpjklm) 

+ {PjkPilmn + PjlPikmn + PjmPikln + PjnPiklm) + (pklPijmn + PkmPijln + PknPijlm) 
+ {plmPijkn + PlnPijkm) + PmnPijkl + 4 (pijPfcZpmn + PijPkmPln + PijPfcnPim) 
+ (PifcPj/Pmn + PikPjmPln + PifcPjnPhn) + (Pi/Pjfcpmn + PilPjmPkn + Pi/PjnPfcm) 
+ (pimPjfcPin + Pimpjlpkn + PimPjnpkl) + (pinpjfcPZm + pinpjlpkm + PinpjmPkl) 

It follows from Theorem 6.4 later that the general formula is 

p6P P] ^en(p) 

where PA is the set of partitions of n using only even numbers and £ the length of a given 
partition p. As usual, H(p) is the set of all disjoint partitions of the set of indices into subsets 
of lengths given p. We denote the resulting subsets of ir by m , . . . , iri and let p^ be the product 
of p- functions with indices the subsets 7Tj. 
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6 7£-functions 

In this section we study a new class of functions, designed to generalise the p and Q-functions 
in turn. Hence we have chosen to call them "/^-functions. 

Definition 6.1. We define the n-index rath order IZ-functions from the Kleinian er-function, 
using the generalised Baker-Hirota derivatives of Definition 3.2, as 



<U = (-z^) Son^o...ont ] &A, (6.i) 




n, -i«n ., no- 
where i\ < ■ ■ ■ < i n € {1, . . . ,g}. 

Like the Q-functions, the indices of 7£-functions do not indicate differentiation, but rather 
the application of operators. 

Corollary 3.5 shows us that the n-index 7£-functions are zero unless n is a multiple of m. We 
can also observe from Corollary 3.6 that for such n the functions are not identically zero. From 
the definition it is clear that the mth order 7£-functions have poles of order at most m, occurring 
on the B-divisor. The main result of this section will be a proof that these functions are also 
Abelian, and hence can be used to complete the basis for T(m). Following the approach taken 
for the Q-functions in the last section, we will do this using the Leibniz rule for the operators 
and the effect of the operators applied to h, the extra factor introduced in the quasi-periodicity 
condition (4.4) for the u-function. 

Lemma 6.2. We have, for all n > 1, 

m 

fe=i 

Proof. We note that if n is not a multiple of m then the result follows directly from Corol- 
lary 3.5. In Lemma 5.2 we could show directly that the symmetrization of the original Baker- 
Hirota operators evaluated on h gave zero due to well known properties of the binomial coefficient 
involved. We saw in Section 3.2 that for the general operators the coefficients involve not just 
multinomial coefficients but also Mobius functions and so we use a different approach. 
First we recall Lemma 3.3, which allows us to conclude that 

m / m \ 

H^o...ont ] ^h = Yjc in (^h^). (6.2) 

Here £ n is an index and the sum is finite with the number of terms depending on n. Each ix\~i 
is a subset of the indices I n = {ii, ■ ■ ■ , i n } such that iti/, ■ ■ ■ , 7Tk,e are a disjoint partition 
of I n . We further recall equation (5.4) which stated that for the function h we have lijj ^ = 
Xe L Li 1 • • • Lj n . Hence each of these tensor products will symmetrize to give the same function, 

m mL t ... r,. 
l / fe=i 

So we have 





S o yW o • • • o ^ h = X m e mL L tl ■ ■ ■ L in [ £ c in | . m , 

fe=i 
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We aim to show that the sum of constants Q n in equation (6.3) is always zero, which would 
allow us to conclude the lemma. For the case n = 1 we have 

mm m 

k=l j=l k=l 

and so l\ runs from 1 to m with c^ = C . Clearly the sum of these constants is zero since C 
is an ?Tith root of unity. 

Now consider applying another Hirota operator to equation (6.2): 

mm mm m 

< ] o-o «£ <g> ft = £ <r ^ E ^ ® k. = E c'- 1 E c ^ <8> ^ • 

fc=i i=i £„ fe=i i=i £ n fc=i 

Since 

""7Tfc>*n+l = ^in + l^k 

we have that for all j = 1, . . . , m, 

(m \ 

5 ?J +1 (2)V.j=X m e mL L 4l ...L Wl . 

So we see that 

m / m 

5 o «H o • • • o ^ (g) fc = x ™e^L n • • • L,„ +1 £ C^'- 1 E ^ 

fc=l \j' = l ^n 

Hence 



7ft 



E^ +1 =Ec^E^ = o- ■ 

We can now prove the periodicity of the "/^-functions. 

Theorem 6.3. The n-index rath order IZ-functions are all Abelian. 

Proof. By Definition 4.6 of the lattice operator we see that the Abelian property is equivalent 
to showing 

A e oTZ lm] ■ =lZ [m] ■ . 

Using the definition of the 1Z- functions in equation (6.1), Lemma 4.7 which concluded that Ap 
commutes with the other operators, and the periodicity property of a, we have 

A e o K l ™ ] % = ( - 1 ^ S o uY" o • • • o nY oA e (fta, (6.4) 

1 «i,...,i n i mh m o~ m I "O v ' 

v J k=\ 

where h = \e . We observe that under the multiplication operation for tensor products 

(m \ m / m \ / m \ 

(g)<r =0/iff= (g)/i (g)a . 
fc=i / fe=i \fe=i / \fe=i / 
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So we may use the Leibniz rule from Corollary 3.9 to give 

mm n / m \ / m \ 

*ir I o--o7if I ®*® - = E£*S ,] W* ®a . (6.5) 



fc=i fe=i £=o Tren \fc=i / \fe=i 



When £ = we see there is only one ir £ II and that 7Ti is the empty set. So the inner sum in 
equation (6.5) produces just one term, 



in- 



\fc=i / \fc=i / 

For all ^ > the inner sum in equation (5.6) produces terms which all involve £ operators applied 
to <8>™ =l h. So from Lemma 6.2, only the term from £ = will not vanish under symmetrization. 
Hence we have 

(m \ / m \ / m \ 

<g>* ®* =no«Ho...o H H (g) a 
fe=i / Vfc=i / \fc=i / 

and so equation (6.4) becomes 



/ 1 \ m 

A e oK [ ™ ] t =( -± U m 5o^ n o...o^ in (S?) C T = i?i ml 

v y fc=i 



Since the "/^-functions are Abelian we know they can be written algebraically using the 
Kleinian p-function. We make this relationship explicit in the next theorem. 

Theorem 6.4. The n-index mth order IZ-functions may be expressed in Kleinian p-functions 

as 

<U = E (-™)'" 1 E **> ( 6 - 6 ) 

where Pn is the set of partitions of n using only integers divisible by m and £ the length of 
a given partition p. As usual, H(p) is the set of all disjoint partitions of the set of indices into 
subsets of lengths given p. We denote the resulting subsets of it by tvx,. . . ,tt£ and let p n be the 
product of p-functions with indices the subsets 7Tj. 

For example, in the m = 3 case we have 

+ Pii,ii,i&Pi3,i4,,i5 ' Pii,i3MPh,i5,i6 ' Pii,i3,i5Pi2M,ie ' Pii,i3,i6pi2M,h 
+ Pii,i4,i5pi2,i3,ie ' PiiM,iePi2,i3,is ' Pii,i5,i6pi2,i3MJ • 

Note that equation (6.6) clearly reduces to equation (5.7) in the m = 2 case. 

Proof. In this proof we define a symbol p such that a = e p . Then p = log(<r), pi = Gi/a and 

Ph,...,in = ~Ph, — ,i n 

for n > 2. It is important to note that p is not the Weierstrass p- function and neither p nor pi 
are Abelian functions. We consider an mth order %-operator acting on a tensor product of a. 

m m m m 

«N a = rl [ : r ; ] (g) e* = £ C'^df <g) et 

fc=l fc=l j=l fc=l 
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This will be a sum of terms, each a tensor product in which every entry has a factor e p . So 
using the multiplication for tensor products, we have 

m I m \ / m \ 

fc=i y=i / \fc=i / 

with the <pi x occurring in the jth position. 
For any subset I & I n = {i±, . . . , i n } define 

m 

£/ = £ C |/|(j_1) (1 ® 1 ® • • • <8) pj <8> 1 <8) • • • <8> 1) ■ 

i=i 

Then equation (6.7) becomes 

(m \ / m \ 

®H=E { « l} (®aj- (6.8) 

Now consider applying a second operator. Using the product rule for %-operators from Lem- 
ma 3.8 we have 

(m \ / m \ / in \ 

®°) =^r ] (s M ) ®* +% } ?4r 1 ®a . (6.9) 

fc=i / \fc=i / \fc=i / 

Applying the ^-operator on Er^i will produce a sum of terms, many of which will have zero as 
one of the entries in the tensor product. In this proof we will be applying operators and then 
symmetrizing, meaning that these terms will all vanish in the end and so we need not keep track 
of them. We let the symbol = indicate not a true equality, but that the difference consists only 
of terms with zero as one of the entries in the tensor product. When we apply the ^-operator 
on S/ the only terms which will not vanish are those where the differentiation is applying to 
the entry with pj. We hence have 

^: ] (s/)=s {/A} . 

This result, along with (6.8), simplify equation (6.9) to 

/ m \ / m \ / m \ 

[ml n i [ml 



\fc=i / \fe=i / \fc=i / 

An induction gives 



771 



«No...o^J ®,=E £ XV <g)a , (6.10) 

Vfc=i / peP„ 7ren(p) \fe=i / 

where Pn is the set of partitions of n and the other notation as stated in the theorem. 
Now consider the application of the symmetrization operator. We have 

(m \ I m 

®o\=a m and 5(E/) = 0/ £ C |J|y_1) ) 
fc=l / \j=i 

As in the proof of Lemma B.3, if m divides \I\ then we obtain 

S(£/) = mpi = -mpi. 
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and obtain zero otherwise. Hence (6.10) becomes 

So*Ho...o«H((g)a)= £ ( _ m) < £ ^ 
Vfe=i / pgP H ^(p) 

and dividing by —ma m gives the result of the theorem. ■ 

It is clear from their definitions that the Q-functions are the 2nd order 7£-functions: 

Qh,i2,...,i n (u) = nfl i2 ^ in {u). 

Hence in turn the two-index p-functions are the 2-index 2nd order 7£-functions: 

Pij(u) =Tlfj{u). 

However, we see from the following corollary that in fact all the Kleinian p-functions are them- 
selves 7£-functions. 

Corollary 6.5. The m-index mth order IZ-functions are exactly the m-index p-functions, 

Pii,ia ) ...,t ra (ti)=^! a> .„ jira (ti). 

Proof. This is the application of Theorem 6.4 in the case n = m. ■ 

Recall that the n-index p-functions had definite parity indicated by n; i.e. they are odd if n is 
odd and even if n is even. This property extends to the 1Z- functions, and hence the Q-functions 
also. 

Lemma 6.6. The n-index 1Z- functions have parity equal to that ofn, with respect to the change 
of variables u — > [— l]u. 

The proof is simple but lengthy by considering the different combinations of parity of the 
components. The full details are in Appendix C. 

7 New bases of Abelian functions 

In this section we present some new bases of Abelian functions, completed using 7£-functions. 
These are for the two canonical genus three cases. In this section we work with functions of 
only three variables and so we drop the commas between indices for brevity. Since the genus is 
fixed so are the dimensions of corresponding bases. However the actual bases themselves require 
different functions, made apparent by the different weights in the two cases (see Section 4.2). 
We note that like the p-functions, the weights of the 7£- functions are the sum of the weights of 
the variables indicated by the indices. This can be seen by considering 7£-functions as a sum 
of terms, each with a m in the denominator and a monomial of cr-derivatives in the numerator, 
which includes all the indices. The weight of a-derivatives is the weight of a minus the weight 
of the indices and in each term the weight of a cancels leaving the overall weight determined 
only by the indices. 

For each T(m) a basis is constructed systematically as follows. We start by including the 
basis for the preceding case so we then need only search for functions with poles of order m. We 
consider a set of possible basis functions, say 7£-functions with a given number of indices, and 
proceed in ascending weight order. We use the a-expansion for the corresponding (n, s)-curve 
to see if a given function is linearly independent of the existing basis entries and if so add it 
to the basis. It is important to include the functions from the preceding bases when testing 
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linear independence, since although the mth order ^-functions have poles of order at most m, 
they may in fact have lower order poles when associated with certain curves and this can only 
be checked using the expansion. Once the set has been exhausted, we consider another set if 
necessary. 

We use the weight theory to consider only linear combinations of functions (and curve pa- 
rameters) which give the correct weight. We then substitute in sufficient terms of the expansion 
and see if there are solutions for the coefficients. Such calculations can be lengthly and grow in 
CPU time and memory requirement with both the genus and the number of poles. Significant 
computational simplifications can be made by writing procedures that take further advantage of 
the weight structure present in the theory. For example, when expanding the product of series 
it is only necessary to multiply those terms which will give the correct final weight, as the other 
terms must all cancel. We can further reduce calculations by noting that there is a finite weight 
range that basis entries can take. Entries in T(m) can have weight no lower than — mwt(cr) (see 
Lemma 3.4 in [16]). 

The bases presented below are all for functions associated with general (n, s)-curves as in 
equation (4.1). However, we only need use the series expansions associated to the corresponding 
cyclic curves in equation (4.2) to check linear independence. If an element cannot be expressed 
using the basis with the restriction on the parameters, then neither will it be expressible with 
the wider set of parameters. Further, we only need to use sufficient expansion to give non-zero 
evaluations of the functions considered in order to check whether they are linearly independent. 

We note that while these bases are not unique, (they can formulated using different functions), 
their weight structures are. Also, there are other invariants of the bases, such as the number of 
odd and even functions. We have several conjectures on these invariants, which are still under 
investigation. 

7.1 Hyperelliptic curves 

The hyperelliptic curve of genus 2 is a special case as its bases can be finitely generated by 
differentiation (see Table 1). If desired, they could be reformulated using 7£-functions. All the 
p- functions are themselves 7£-functions and the extra function, A £ T(3), may be replaced by 
the six index 7£- function, ^122222- Of course > A may be expressed as a linear combination of 
the basis entries. In the cyclic case the relationship is, 

A- 2 \ -7?[ 2 U 4 \ 2 7?[ 2 ] 7 \ 7? PI 1 V^ 9 \ 

~~ _ 5 u 5 12 ~~ 5 12 ~ 20 122222 ~ ^ A i- 

In the genus 3 hyperelliptic case we find that a basis for T(2) is 

(1 K [2] K [2] TZ [2] K [2] K [2] K [2] K [2] \ (7 V) 

X 1 : /v ll> /v 12' /v 13' /v 22' /v 23' /v 33' /V 2222J • \' ■ ± ) 

f2l 

We note that 7^ 2 2 2 2 ^ s a second order ^-function, and so equal to Q2222 as defined in Section 
5. We note further that there are two other 4-index Q-functions of this weight; Q1133 and Q1223 
and either of these may be used instead. Indeed, there are even more alternative Q-functions 
that could play the role if we allow more indices. 

We now proceed to give a basis for T(3). We arrange the elements in the set as vertical lists 
purely to make the structure (number of functions of each type) clearer to the reader 



(7 -1 \ -7-. [3] 7-. [3] v [3] v [3] v [3] o 7-, [2] 

{<■■>-), 'M13' 'M33> 71 -233' '^233333' '^222222' c, l /t 2222' 

^[3] ^[3] ^[3] ^[3] -,[3] ^[3] rl1 -,[2] 

'Mil' / ' c 122' / ' c 222> /,c 333' /-t "133333' 'M 13333' c, 2'v 2 222> 

^[3] ^[3] ^[3] ^[3] ^[3] fl-pp] 

"'112! / ' c 123' / ' c 223' /-t "222333' / ' c 112333» c '3'^2222 



(7.2) 
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Here (7.1) is representing the eight entries in the basis for T(2) presented in equation (7.1). We 

[21 

note that the three first derivatives of ^2222 cou ld be replaced by 9- index third order 1Z- functions 
if desired. Finally, a basis for T(4) is 



(7-2), 
'Mill' 

'Mll2' 
'Mll3' 
'Ml22' 
'Ml23' 
'Ml33' 
'M222' 



'M223' 
'M233' 
'M333' 
'M222' 
'M223' 
'M233' 
'^2333' 
'^3333' 



d 3 n 

d 2 K 



[3] 

233333' 

[3] 

233333' 
a -r> [3] 
^3 't-133333 ) 

a t?[3] 

"1 / ' c 233333 ' 

^2^133333, 

fj -t?[3] 

u l /-t -133333 ' 

fj -t?[3] 

a 2 '^222333' 

fi -r> [3] 

"1 / ' c 222333 ' 



fl -jp [3] 
^'M 13333' 
a -i?[3] 
^l'M 13333' 

Pi ipPl 
U 2'M. 12333' 



»[3] 
"222222' 

^1^-112333' 

,[3] 

"222222' 
a -15 [3] 
c, l' v 222222' 



W 



<9 2 ^ 



'M2333333' 
'M1333333' 
'M1233333' 
'M2223333' 
'M1133333' 
'^22222222' 
'M1123333 



> . 



In this case those first derivatives of third order 7£- functions could not be replaced by any fourth 
order ^-functions. We may conclude this because previous calculations in [16] have shown that 
the basis must contain some odd functions with poles of order 4. However, since the number of 
indices on the 4th order 1Z- functions must be a multiple of 4, all these functions must be even 
by Lemma 6.6. The odd functions with poles of order 4 are all found within the derivatives of 
the previous basis. 

Although we have not proved that all bases in this case can be constructed this way, we have 
overcome one of the key difficulties in constructing such bases, which becomes apparent when 
examining the weights. In this case the variables u = (ui, 112,113) have weights (5,3,1) while 
the a- function has weight 6. A basis for T(m) must have an entry with the minimal weight of 
— mwt(cr) (see Lemma 3.4 in [16]). So in this case each T(m) requires a function of weight —6m 
in its basis and such a function cannot be a derivative of a function in T(m — 1). This is since 
such a function could have weight no lower than 

— 6(m — 1) — 5 = —6m + 1 

achieved by differentiating a minimal weight function in T(m — 1) by U\. We can now always 
choose the minimal weight function in the basis of T(m) to be the 2m- index function IZ^ 2 - 



7.2 Trigonal curves 

The trigonal curves are those (n, s)-curves with n = 3. The simplest is the (3, 4)-curve which 
also has genus three, but is non-hyperelliptic. We find that a basis for T(2) is 



r, v [2] v [2] ^[2] v [2] v [2] ^{2} ^ [2] 1 
X 1 , 'Ml ' /v 12' /v 13' /v 22' /v 23' /v 33' ' V 2222J- 



(7.3) 



We note that symbolically, this looks the same as equation (7.1). However, the two theories 
differ in the weights of the functions. In the hyperelliptic case u = (111,112,113) had weight 
(5,3,1) while in the trigonal case the weights are (5, 2, 1). 

[21 

Once again, 7£ 2 222 = Q2222 and there were other choices of Q-function. This time there was 
only one alternative 4-index Q-function of this weight, Q1333, contrary to the alternatives in the 
hyperelliptic case, making the symbolic equivalence somewhat coincidental. However, there is 
also a structural equivalence in the use of six 2-index functions and one 4-index function. 

For r(3) in the trigonal case we find the basis 



(7.3), 


TZ [3] 

'Ml3' 


TZ [3] 

'M33> 


'^233' 


'^333333' 


'Ml3333> 


Ft "7? PI 

a l /v 2222> 


'Mil' 


'M22' 


/v 222' 


'^333' 


'*"222333> 


TZ [3] 

'M22333' 


<^2'^-2222' 


'Ml2' 


'M23' 


TZ [3] 

'M23> 




K [3] 

'M33333' 


'M22233' 


f) -7?[ 2 1 
<y 3'M222> 



> . 



(7.4) 
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So the symbolic equivalence is now lost; for example the function ^333333 is present in this basis 
but can not replace any function in the basis (7.2). However, the structural equivalence is still 
present; i.e. in both the hyperelliptic and trigonal cases the number of functions with a given 
number of indices was the same. Similarly, for T(4) in the trigonal case we find the basis 



(7-4), 
'Vlill> 

/v 1112' 

'^1113' 
'^1122' 
'M1231 
'Ml33' 
/t 1222! 



'M223' 

TZ [A] 

'M233' 

'M333' 
/t 2222> 

'^2223' 
'^2233' 
'^2333' 
'^-3333' 



d 3 n [ 3 
o 2 nl 



.[3] 

"333333' 
,[3] 
"333333' 

^3^222333' 

fj -t?[3] 

a l /-t -333333 ' 

d 2 n l ? ] 



d 2 n\ 



'222333' 

,[3] 

'133333' 

J3] 

-122333' 

^1^-222333' 



d 2 n\ 
d 3 n [ ; 



,[3] 

'122333' 
a t?[3] 
a l /,c 133333 ' 

^2^122233' 

a -T?[3] 
C '2'M 1 3333, 

f) -r> [3] 

a l /-t "122333 ' 

f) -t?[3] 

a l /-t "122233 ' 

f) -r> [3] 

'-'I'Missss' 



ft [4] 

'^22333333' 
'^-22233333' 
'^13333333 > 

'M2233333' 
'M2223333' 
'^11333333' 
'^11223333 



> . 



Again, the number of functions of each type is the same. We note that although bases for these 
spaces have been derived before in [18] and [16], this structural similarity only becomes apparent 
now, when using 7^-functions. 

We conjecture that in general a basis for T(m) can always be constructed using the preceding 
basis, rath order ^-functions and the first derivatives of (m — l)th order 7£- functions, which 
would solve the problem of identifying enough suitable basis functions. 

A Proof of Lemma 2.3 

We need to verify equation (2.1), 

VhO---oV in {f®g) = ]T(-ir E /M^.%...-,^}®^!^.-^ 

m=0 y l<ji<J2<— <jm<n 

Proof. We give an inductive proof. First we observe the case n = 1 where 



RHS(2.1) = (-1)° (/ M0 ® g ) + (-1) 1 Y, fh\H h } ® 9i h 

= fh ®90-f0®gii = fh®g-f® 9h 

as required. We now assume the result holds for n = k and consider the case n = k + 1. 
V h o ■ ■ ■ o V ik+1 (f®g)= V ik+1 o V h ■ ■ • o V ik (f g) 

\m=0 \l<ji<J2<---<jm<A; 

m=0 y \l<ji<J2<---<jm<fc 

2_^ Jl k \{i jx ,i^,...,ij m } ^ 9i jl ,ij 2 ,...,i jm ,i k +l 



(A.l) 
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We need to show that equation (A.l) leads to 

V h o---oV ik+ ^f®g) = J2(-irl Y, fh + i\{in,i^^rn^9i h ^ 2 ,..,i jm ) -(A.2) 

m=0 \l<jl<J2<--<j m <k+l J 

Each value of m in the outer sum in equation (A.l) gives two parts. The first part will contribute 
to one value of m in equation (A.2) with the second contributing to the following value. These 
two parts are required to fill each value of m in equation (A.2), with the exception of the first 
and last terms in the sum, for which one part is sufficient. 

Consider the contribution of the first part of equation (A.l) when m = 0. We obtain fi k+1 ®g 
which is what we require for m = in equation (A.2). 

Similarly consider the contribution of the second part of equation (A.l) when m = k. We 
obtain (— l) fc+1 / (g) gi k+1 which is what we require for m = k + 1 in equation (A.2). 

Now we consider the general term. We combine the contribution of the first part of equa- 
tion (A.l) when m = m with the contribution of the second part when m = m — 1 to obtain 

(-!) m I 2^ fl k+1 \{i h ,i h ,...,i jm } ® 9i h ,i J2 ,...,ij m 

\l<jl<J2<-<jrn.<k 

+ Z^ fh\{i] 1 ,i] 2 ,-,ij m _ l }®9i Jl ,ij 2 ,...,i Jm _ 1 ,i k+l 

i<ii<i2<--<im<fc 

Note that the minus sign in the second part was combined with the (— l) m_1 outside to factor 
a (— l) m out altogether. 

We can check that the terms here are referring to the same partition of the set of indices over 
the tensor product. Counting in the first part we see we have k + 1 — m indices and m indices. 
Similarly in the second part we have k — (m — 1) indices and m — 1 + 1 indices, the same values. 

Now note that the terms in the second part never contain ik+i on the left of the tensor 
product, but always contain it on the right. Hence these terms are those extra ones we would 
obtain if we ran the sum in the first part to k + 1 instead of A;. So these terms can be combined 
into the one sum to give, 

(-!) m I 2^ fi k+1 \{i J1 ,i J2 ,...,i Jm }®9i Jl ,i J2 ,...,i j 

\l<7'l<J2<-<jm<fc+l 

the term in equation (A.2) when m = m. ■ 

B Symmetric functions 

We gather together some definitions and results for symmetric functions which we require. 

Definition B.l. The monomial symmetric polynomials in m variables, x = [x±, . . . ,x m ] are 
each defined by a sequence of natural numbers of length m, say p = [pi, . . . , p m ], as 

M p (x)= Y, x ^- 
V>e*(p) 

Following the notation introduced for Lemma 3.3, ^(p) is the set of all permutations of p with 
entries denoted by p = [px, . . . , p m ]. We then let x^ denote a corresponding monomial, 

~ll> — _^l „>2 . . . Vm 
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We define the augmented monomial symmetric polynomials as 
M p (x) = Y, ^ P - 

This time the sum is over permutations of the variables rather than permutations of the powers. 
For example, 

M[3,1,1](X1,X 2 ,X 3 ) = x\ X 2 X 3 + X1X2X3 +X1X2X3, 

^[3,1,1] ( x li x 2i X s) = xfx 2 Xs + xfx S X2 + x\x\X 3 + X^X^Xi + x\xiX 2 + x\x2X\ 

= 2M[ 3il>1 ](xi,X2,x 3 ). 
Let r]i be the multiplicity of % E p, i.e. the number of entries in p equal to i. Then we have 

M p (x)= (i\ V Am p (x). (B.l) 

The monomial symmetric functions naturally occur when considering the application of H- 
operators in Section 3.2. They occur evaluated on the roots of unity. Symmetric polynomials 
in the roots of unity have been considered previously in [30] , with nice results for some of the 
other classes. For example, the elementary symmetric polynomials reduce to 

ei = 0, e 2 = 0, . . . , e m _i = 0, e m = (-l) m+1 . 

However, in our case, the power sum basis is the most useful. 

Definition B.2. The Newton power sums in m variables xi, . . . , x m are defined for k > as 

Vk\X\, • • • , x m ) = Xi + x 2 + • • • + x m . 

Any symmetric polynomial in m variables can be expressed as a polynomial with rational 
coefficients in p±, . . . ,p m . We observe how they evaluate on our choice of variables. 

Lemma B.3. The power sums evaluated on X = [X\, . . . , X m ] where X^ = ( k ~ 1 for ( a primi- 
tive mth root of unity satisfy 

{m ifrn\k, 
otherwise. 

Proof. We have 

1 /-km 

Pk (x) == 1 + c k + C 2k + ■ ■ ■ + ( {m - 1)k = -ttV- 

Since £ is an mth root of unity, the numerator is always zero. The denominator is only non-zero 
for m\k. Hence if k is not a multiple of m then pk=0. If k = m then we have 

p k = X™ + X 2 m + • • • + X™ = l m + C m + C 2m + • • • + C (m_1)m = m. 

Similarly if m\k then we have a sum of m terms, each £ raised to a power divisible by m. Hence 
the polynomial evaluates to m. ■ 
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We now introduce a procedure for expressing symmetric functions in power sums. Although 
every symmetric function of m variables may be expressed as a polynomial in the first m power 
sums pk, this explicit procedure will instead give an expression in the first n power sums, where n 
is the number partitioned by the sequence p = [pi, ... , p m ] . 

We again let II be a set of partitions of indices. This time it will be all disjoint partitions, n, 
of the set of indices I m = {i±, . . . , i m } into subsets denoted by tt\, ... ,n£. Note that we do not 
restrict to a given number of subsets or specify their size. We see that £ will range from 1, 
(when 7r has just tt\ = I m itself), to m (when it consists of the m subsets with one index each). 
We define the associated Mobius function ^(tt) by 



MvOH-ir^ntw- 1 )' 



j=l 



and refer readers to [35] for more information on such functions. 
Next we define length £ vectors v = \v\,...,vi\ 6 N e by 



u(p,n) 



2 Pi'- 



> Yl pi 



Finally we define the function p u as the product of power sum functions whose indices are 
determined by the values in the vector v, 

Pu = Pv\ ' ' " Pu2 ■ 

We can now state the following formula of Doubilet from [15], 

M P = 5^/*(t)Pi/(p iW )- 

wen 



(B.2) 



So the augmented monomial symmetric functions have an integer polynomial expression in the 
power sums p\ , . . . , p n . The expression for the monomial symmetric functions will be the same 
but with a fractional multiplicative constant depending on the multiplicities of p as described 
in equation (B.l). 

We give an example in the case m = 3. We consider the disjoint partitions of ^3 = {i±, 12, is}- 
The set II has five elements: 



7T 



[1] 



TT 



TT 



TT 



TT 



TT 



[1] 



TT 



TT 



TT 



TT 



TT. 



TT. 



TT. 



TT. 



M, 



TT. 



[5] 



{is}. 



{h,i2,h}, 

{h,i2}, 

{h,h}, 

{i2,«3J, 
'1 -{*l} 

The Mobius function then gives 

/U (ttW) = (-1) 3 - 1 2! = +2, 

li(n [2] ) = p{tt [3] ) = p(ttW) = (-1) 3 ~ 2 1!0! = -1, 

/U ( 7r [ 5 ]) = (-l) 3 - 3 (0!) 3 = +l. 

Hence any augmented monomial symmetric function in 3 variables may be expressed using the 
power sums in 3 variables as 

Mp = 2pp 1+ p 2+P3 — pp 1+P2 p P3 — pp 1+P3 p P2 — p P2+P3 p pl — p Pl p P2 Pp 3 - 
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For example, 

^[3,1,1] = 2 P5 - 2p 4 Pl - P2P3 +P3P1- 

Note that when we have entries in p which are zero then we need to recall that po = m. So in 
the m = 3 case we have, for example, 

A%o,o] = 2 P3 - 3p 3 Po +PsPo = 2p 3 - 

C Proof of Lemma 6.6 

We prove that the ^-functions have parity equal to that of n, the number of indices. 
Proof. Think of the n-index mth order 7£-functions as 

r m ] polynomial in u-derivatives m t\ 

K ii,...,i n = ~^ • V-" 1 ) 

Each term in the polynomial in the numerator is a monomial in a and its derivatives containing 
re-indices overall (see Lemmas 3.3 and 3.4). Recall that the cr-function is either odd or even. 
In the case of (re, s)-curves this can be determined from (n,s) (see Lemma 4.4). We refer to 
derivatives of a with respect to an odd number of variables as cx -derivatives and similarly 
<r e -derivatives if differentiated with respect to an even number of variables. 

Suppose first that a is an even function and so the denominator of (C.l) is even. In this 
case the <r -derivatives are odd and the cr e -derivatives even. If n is odd then each term in the 
numerator of (C.l) must contain an odd number of <r -derivatives. Similarly if n is even then 
each term would have an even number. Hence the numerator, and therefore 7Zi „• , has parity 
the same as n. 

Now suppose that a is an odd function. This time the denominator of (C.l) has parity equal 
to that of m, and it is the <7 e -derivatives which are odd while the <r -derivatives are even. In this 
case we need to be more careful. Note that each term in the numerator of (C.l) can be thought 
of as a product of cr-derivatives multiplied by the cr-function itself raised to a power. 

If n is odd then m must be also, making the denominator of (C.l) odd. Each term in the 
numerator must contain an even number of a e - derivatives, paired with an odd number of a Q - 
derivatives. Hence, since m is odd, any cr-functions in the term must be raised to an even power. 
So overall each term in the numerator is even and hence 1Z™ 1 „■ is odd overall. 

If n is even then we must consider the two choices for m separately. If m is even then so 
is the denominator of (C.l). Each term in the numerator must have an even number of a- 
derivatives overall. If there are an odd number of cr e -derivatives then there are an even number 
of <T -derivatives and hence, since m is even, a cr-function raised to an odd power. Then each 
term in the numerator is even overall. If there were an even number of <r e -derivatives then there 
are an odd number of do-derivatives and any remaining cr-function is raised to an even power. 
Either way, each term in the numerator is even overall and hence so is 1Z™ i . Finally, if m 
were odd then the same analysis will this time show that each term in the numerator is odd. 
However, the denominator will also be odd in this case and so TZ;- ■ is still even overall. ■ 
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